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ABSTRACT 

Single  component  plasmas  are  found  to  exhibit  a  reversal  In  the  roles  of 

cyclotron  and  upper  hybrid  resonances  familiar  In  neutral  plasmas.  There  Is  a 

'\  ' 

global  resonance  at  the  cyclotron  frequency  u  >  0  =  qB/mc  and  a  local  field 


cancellation  at  layers  where  u  *  Qj  »  ojp  being  the  local  plasma 


frequency.  Thermal  effects  limit  the  amplitude  of  the  global  resonance  and 
produce  energy  absorption  at  the  u  «  0]^  layers. 


Aoor’sslon  For^ 
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PACS  numbers:  52.25Wz,  S2.20Dq,  52.35Fp 


single  component  plasmas  confined  by  external  magnetic  fields  strong 
enough  to  overcome  the  Coulomb  self-repulsion  have  been  demonstrated^*^  to 
exhibit  remarkably  long  lifetimes  and  to  be  rather  quiescent  compared  to  other 
laboratory  plasmas.  Such  properties  make  these  systems  ideally  suited  for  the 
study  of  effects  in  which  the  presence  of  a  self-consistent  electrostatic 
field  plays  an  important  role.  A  proper  understanding  of  wave  and  transport 
phenomena  Including  the  effects  of  such  self-electric  fields  is  essential  for 
developing  a  better  description  of  conventional  "neutral"  plasmas,  which  in 
reality  develop  strong  amblpolar  electric  fields,  due  to  the  different 
transport  properties  of  electrons  and  ions.  The  present  theoretical  study 
examines  some  unusual  properties  of  an  extreme  nonneutral  system,  namely  a 
single  component  plasma,  with  special  emphasis  on  the  plasma  response  near  the 
cyclotron  frequency.  The  effects  analyzed  may  also  be  of  Interest  in 
understanding  cyclotron  resonance  absorption  and/or  emission  in  plasmas 
confined  in  magnetic  mirror  configurations,  where  electrostatic  fields  are 
known  to  be  present. 

In  a  cold  magnetized  neutral  plasma,  cross-field  plasma  behavior  near  the 
electron  cyclotron  frequency  ^  is  governed  by  the  dielectric  tensor  element 

-  1  -  ai^/((D^-n^)  •  ((ii^-aij^j^)/(a)^-n^)  where  Up  is  the  electron  plasma 
frequency  and  +  (0^)  the  upper  hybrid  frequency.  Two  features 

of  e^x  chat  have  important  consequences  for  the  plasma  response  are 
e^x  *  ^  3is  u)  n  and  Cxx  0  as  o)  For  a  neutral  plasma  slab  in  which 

the  gradient  in  the  density  OqCx)  is  perpendicular  to  a  uniform  magnetic  field 
Bz,  the  application  of  a  capacitor  plate  field  E^x  at  a  frequency  to  results  in 
a  collective  field  which  vanishes  at  all  x  when  oi  ■  .T  and  has  a 

local  resonance  (Ej  +  *)  at  the  upper  hybrid  layers  where  m  «  a)yjj(x). 


For  a  single  component  plasma  we  find  a  strikingly  different  behavior. 

The  dielectric  tensor  element  Cxx  ^  cold  single  component  plasma  can 
be  obtained  as  the  lowest  order  term  in  a  kinetic  treatment  which  Involves 
integrating  the  linearized  Vlasov  equation  along  the  orbit  of  a  particle 
moving  under  the  Influence  of  the  self-electric  field  Eq(x)x  and  the  uniform 
magnetic  field  Bz.  Assuming  the  cyclotron  radius  r^  «  L  =  l(d/dx)ln  no 
and  retaining  only  terms  linear  in  to/I*  in  the  electric  field 

Eo(x')  =  EqCx)  +  dEo/dx  (x'-x)  »  EqCx)  +  4iTqno(x)(x’-x)  (1) 

governing  the  orbit  x'(t)  of  a  particle  starting  from  x,  we  find  that  the 

linear  term  in  the  expansion  tends  to  weaken  the  magnetic  restoring  force 
and  the  particle  gyrates  not  at  angular  frequency  but  at  H^Cx)  = 

[Q^-(ij^(x)  Hence,  for  a  single  component  plasma,  Ejqc  =  1  -  =• 

Comparison  with  the  neutral  plasma  result  shows  that  the 

familiar  roles  of  single  particle  resonance  and  collective  resonance  are  in  a 
sense  Interchanged;  we  now  have  Cjm  ♦  0  as  u  JJ  and  c^x  “  at  layers  where 

oa  =  n^Cx)  »  ^ The  collective  field  Ej  exhibits  a  global  resonance 

(l.e.,  El"*-*  for  all  x)  when  oi  ■  n  and  a  local  null  at  layers  where 

ti)  »  J2  (x). 

1 

As  expected,  finite  temperature  effects  modify  both  the  global  resonance 
and  the  local  null  of  Ep  We  now  discuss  two  such  effects.  The  first  one, 
characteristic  of  nonneutral  plasmas,  arises  from  the  modification  of  the 
single  particle  orbits  when  higher  order  terms  in  r^/L  are  included  in  the 
expansion  (1)  for  the  equilibrium  electric  field  EqCx').  The  second  one  is 
associated  with  the  existence  of  collective  modes  and  can  be  regarded  as  a 


correction  to  e^x  in  powers  of  kr{.  where  k  =  |d/dx  lnEi|,  Ej^  being  the  linear 
wave  field.  To  evaluate  both  of  these  effects,  we  need  to  choose  a  suitable 
equilibrium  distribution  function  fo(*»v)  constructed,  of  course,  from  the 

constants  of  the  motion.  We  find  that  the  choice  fo“C  exp[-(v^+2q!>o(x)/m)/2V^- 
(x  +  Vy/^)^/2i^],  where  7  and  I  are  given  parameters  and  C  is  a  normalization 
constant,  provides  a  convenient  bell-shaped  density  profile,  peaked  at  x«0.  The 
function  fg  is  a  bi-Maxwellian  in  a  local  frame  drifting  with  velocity  vj(x)  » 

which  is  the  sura  of  the  ExB  and  the  diamagnetic  drift  velocities; 

the  mean  squared  velocity  components  are  <Vx>  ®  <Vz>  »  <(vy  -  vj)^>  » 

v^£^/(r^+i^) .  Note  that  the  equilibrium  electric  potential  appearing  in  fg 
must  be  self-conslstently  determined  from  Poisson's  equation  which  can  be  written 

in  the  dimensionless  form  d^'Ji/dC^*  e’^  ~A,  where  ^  *  -q  [i^Q(x)-(tQ(0)  ]/mv‘^-A5^/2 , 

E  =>  x/Xj),  Aq  •  v/Mp(0),  A  *  AjJ/(r^  +  The  parameter  A,  which  determines  the 

density  profile  no(€)/no(0)»€'l',  must  be  greater  than  unity  for  a  bounded  plasma 
[no  +  0  as  ?  +  «].  The  self-consistent  potentials  and  the  corresponding  density 
profiles  for  A«l.l,  2.0  are  displayed  in  Fig.  1.  For  A  in  this  range,  the 
density  scale  length  L  is  of  the  order  of  the  Debye  length  Xj)  and  so  the 
assumption  r^.  <<  L  implies  a)p(O)  «  0;  this,  in  turn,  implies  that  the  single 
particle  resonance  frequency  is  very  close  to  the  cyclotron  frequency  Q 

for  all  X. 

The  first  of  the  finite  temperature  effects  arises  from  the  inclusion  of 
the  (rc/L)^  correction  term  ( l/2)(d^Eo/dx^)(x'-x) ^  in  Eq.(l)  in  calculating  the 
single  particle  orbits.  The  gyration  frequency  of  a  particle  with  initial 
position  X  and  velocity  v  is  0]^,  evaluated  at  the  guiding  center  position 
Xg  s:  X  +  (qEo/m+nvy)/Q'^.  Thus  the  gyration  frequency  is  modified  from 
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ni(x)-[SI^  +  to 

r  2  q  q  rqE-(x)/m+fVy 

'^2(x,Vy)-fn  +^o'(Xg)]  =f?i(x)+ ^o”<x)f - -p - i] 


Solving  Che  linearized  Vlasov  equation  by  integration  over  particle  trajec¬ 


tories,  we  find 


where  Au  ■  r{.l|d(Dp/dx(  [2J5^(rc  and  Z  is  the  plasma  dispersion  func¬ 

tion.^  When  |u  -  Qi|»A(ii,  reduces  to  Che  cold  plasma  result,  but  at  layers 
where  !i)«Qi(x)  it  remains  finite  and  E^  no  longer  vanishes  there.  Furthermore, 
Gxx  is  complex,  corresponding  to  energy  transfer  from  the  field  to  the 
particles.  This  energy  transfer,  which  arises  from  the  velocity  dependence  of 
^2*  is  analogous  Co  Che  familiar  cyclotron  damping  associated  with  the  Doppler 
shift  in  cyclotron  resonance  from  H  to  D  +  when  the  parallel  wave  number 

k2*0.  (In  fact,  allowing  finite  k2  in  Che  present  study  would  result  in 

additional  cyclotron  damping,  which  would  be  negligible  for  k2V  «  Au  or 

k2l>>L(n/(ap)^] .  Fig.  2  displays  the  real  and  imaginary  parts  of  the  electric 

field  E^  -  Ec'^^x  ^*-xx  Ki^en  by  Eq.(3)]  driven  by  an  external  antenna  at 
'o/r!«.9975  for  a  plasma  with  (iip(0)/JJ  ■  .1,  A  ■  1.1.  The  spatial  width  of  the 
resonance  layers  where  Im  E^  is  significant  is  of  the  order  of  r^.,  so  the  power 
absorbed  per  unit  area  of  each  resonance  layer  is  P/a-r^C u/8x)  jEj^l^Im  e^x 

corresponding  to  a  local  absorption  rate  P/ [arcE^/Sir]  ~  (rc/L)n. 


The  inclusion  of  (rj./L)^  corrections  to  the  particle  orbits  removes  not 
only  the  local  null  of  at  layers  where  u  >  but  also  the  global 

singularity  of  as  (except  where  dn^/dx  *  0).  The  global  singularity  in 
E^  at  u)  >  n  is  also  removed  (even  where  dn^/dx  >  0)  by  a  different  thermal 


effect,  namely  the  inclusion  of  finite  kr^.  corrections  to  This  converts 

into  a  differential  operator.  Keeping  the  lowest  order  corrections  whic 

.  ]  L 

Involve  (kr^)^  ,  we  obtain  a  second  order  differential  equation  for  E^, 


where,  in  order  to  emphasize  the  finite  krc  effects,  we  have  suppressed  the  Z 
functions.  Eq.(4)  predicts  the  existence  of  collective  thermal  modes  with  a 

WKB  dispersion  relation  (kr^)^*  (u)^-n^)(a>^4£lj^)/3n^(i^,  propagating  for  oKH, 

(i)>2n]^  with  cutoffs  at  (o^Q,  20^.  These  modes  are  analogous  to  Bernstein  waves 
of  neutral  plasmas.  The  Bernstein  dispersion  relation  for  a  single  species 
plasma  has  been  obtained  before  by  Davidson*^  for  the  special  case  of  a 
cylindrical  rigid  rotor  and  Is  Identical  to  the  neutral  plasma  dispersion 
relation  except  for  the  shift  n  ^  H  -  2(i%  (where  is  the  rigid  rotor  angular 
velocity)  caused  by  the  Coriolis  force.  In  our  slab  problem,  however,  the 
Coriolis  force  is  absent  and  the  shift  Q  -»■  Q|^  Is  caused  by  the  gradient  In  the 
electric  field  (corresponding  to  a  shear  In  the  angular  velocity  for  the 
cylindrical  case). 

The  thermal  mode  electric  field  for  the  undrlven  case  (E^-O)  obtained  by 
numerical  Integration  of  Eq.  (4)  Is  displayed  in  Fig.  3  for  the  choice  ui/fl  « 
.9975,  (Op(0)/fJ».l,  A>>1.1.  (Note  that  the  condition  krj.<l  for  which  Eq.  (4)  is 
valid,  requires  frequencies  close  to  the  cutoff  values.)  In  a  driven  system 
(E(.  *  0),  the  thermal  mode  pattern  appears  superposed  on  the  cold  plasma 
response.  In  a  slab  of  finite  width  (-x^^ 


’(b),  the  thermal  modes  can  be 


strongly  excited  at  those  frequencies  aij,  for  which  the  undriven  solution  has 
nodes  at  x  »  ±X(,.  These  are  the  analogues  of  the  Tonks-Oattner  resonances^  in 
neutral  plasmas. 

Note  also  that  the  existence  of  the  thermal  modes  removes  the  cold  plasma 
global  resonance  at  the  frequency  m  *  Q,  Integrating  Eq.  (4)  twice  for  oi  »  Q, 
we  obtain 


El(x) 


Ec  ^ 


l- 


(5) 


so  that  E^(0)  -  (X{j/r^)^[ai2(0)/a3^(xb)  jEj,  is  large  but  finite. 

In  summary,  the  behavior  of  a  single  component  plasma  for  frequencies  in 
the  neighborhood  of  the  cyclotron  frequency  differs  significantly  from  the 
neutral  case.  New  phenomena  are  available  for  theoretical  and  experimental 


Investigation  that  can  provide  deeper  insight  into  the  nature  of  cyclotron 
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Fig.  1 


Fig.  2 


Fig.  3 


FIGURE  CAPTIONS 


Self-consistent  density  profiles  (top)  and  corresponding  potentials 
(bottom)  in  a  single  component  plasma  for  A  =  1.1,  2.0. 

Real  and  imaginary  parts  of  the  electric  field  Ej  =  Eq/ ^xx  f^xx  given 
by  Eq.  (3)]  driven  by  an  external  field  E,.  for  A  =  1.1,  (i>p(0)/i2  =  0.1, 
u)/f2  =  0.9975. 

Thermal  mode  electric  field  for  the  undriven  case  (E(.  =  0)  for 
A  =  1.1,  u)p(0)/0  =  0.1,  (o/n  =  0.9975. 


Fig.  3.  Thernal  mode  electric  field 
A  -  l.l,  (i»p(0)/Q  -  0.1,  a/a 
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